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Abstract 


ABSTRACT 


In this thesis, an induction motor drive incorporating a fiizzy logic controller to control its 
torque indirectly has been simulated and practically implemented .Tlie present work exploits 
the advantages offered by a fuzzy logic controller (intuitiveness, simplicity, easy 
implementation, and minimal knowledge of system behavior) to obtain a robust, fast, and 
precise closed-loop control of speed of induction motor. Tlie control algoritlim of the drive 
system has been implemented by a pentium PC, which uses PCL-208 data acquisition card 
for communicating with the inverter and the motor. The simulated and practical results show 
that fuzzy logic controller can give fast torque response and precise control of speed of 
induction motor drive system 

Keywords: Indirect field oriented control, Vector control, Fuzzy control. Data acquisition 
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Chapter 1 
Introduction 


1.1 Background 

Alternating cunent iikUdis lu'e getting more and more popularity for applications in industrial 
environment due to characteristics of higher efficiency, liigh torque-to-weight ratio, smatl 
volume and lower cost. The capability to operate at liigher torque, liigher speeds, larger 
power ratings, robust construction and low maintenance cost, make the induction motor 
drives more attractive than dc drives. 

Earlier dc machine had been the obvious choice for applications in adjustable speed drives 
where good torque and speed controllability were desired. The drawback with dc machine is 
that they require ftequent maintenance, have low torque-to-weight ratio, higher volume and 
the sparking at brushes due to mechanical commutation. However, a separately excited dc 
machine is a doubly excited system leading to independent control over both field flux and 
armature current, and so by adjusting these parameters flux and torque can be controlled 
independently. Furthermore, due to low armature inductance the armature cuiTcnt can be 
changed with good dynamic response for fast torque transients. Since induction motor is a 
singly excited system, the flux and torque are not inherently decoupled. So their independent 
control is not possible without a suitable controller. In 1971 Blasclike [1] proposed field 
oriented control or vector control of induction motor by which induction machine can be 
controlled quickly and stably by controlling the stator cun-ent as a vector quantity with a flux 
component and a torque component in two ways. 

(a) Flux vector feedback method where magnetic flux vector from induction motor is fed 
back and stator current is controlled on the basis of flux vector (direct vector control). 

(bj Slip frequency control method [21 where angular velocity of rotor is fed back and is 
added to calculated slip speed to contml the stator current vector (indirect vector 
control). 
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Field oriented control is the most successful in meeting the requirement of wide speed range 
and fast torque response. Using field oriented control, a highly coupled induction motor can 
be controlled tlirough a linear independent decoupled control of torque and flux similar to 
separately excited fully compensated dc motor. 

The basic principle of field-oriented control is to maintain secondary (rotor) flux constant, 
this also assures the control of flux and torque independently. Several schemes exist to 
realize field-oriented control [3-7]. In [7], a torque reference is generated with a calculated 
speed reference and then inverter is switched to track the reference current. The present work 
uses a fuzzy logic controller, wliich sets the torque component of current reference based on 
the speed error and change of speed error. The inverter is switched to follow the set current 
reference witliin a hysteresis band, and this control gives a constant rotor lliix. Fuzzy logic 
controller combines human experience as a rules base so it can account for non-linearity in 
the system behavior into account. Since induction motor also behaves as a non-linear system 
because of change in resistance (due to temperature change) and change in inductaitce (due to 
saturation), a fuzzy logic controller is best suited for controlling an induction motor. 
Furthermore, the fuzzy rules defining the system behavior can be optimized to get best 
possible perfomiance. 

1.2 Objectives 

The present work deals with design, simulation and implementation of an indirect vector 
controlled induction motor using a fuzzy controller. 

1.2 Sununary of chapters 

Tlie first chapter gives a brief introduction and summary of the thesis. The second chapter is 
devoted to computer simulation of the proposed system, which contains the equations 
describing the behavior of induction motor, steps in simulating a fuzzy logic controller and 
algorithm for generation of computer program. The third chapter deals with the realization of 
the proposed system and introduces the hardware used, programming needed to generate 
appropriate signals for control of the drive system. Tie fourth chapter contains the simulation 
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and practical results. The conclusion is drawn in the fifth chapter with a scope for future 
research work. 
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Chapter 2 

Digital Computer Simulation 

2.1 Introduction 

Before implementing a properly designed system, it often becomes necessary to know the 
feasibility of proposed system. In practice, behavior of system is investigated through 
computer simulation. Simulation gives the mathematical solution of equations governing the 
system perfonnance. The solution of performance equations can predict approximate 
behavior of the system. 

In the present chapter, behavior of the proposed system is investigated by mathematical 
solution of equations governing system perfonnance. In the end, some simulated results have 
been shown. In addition, effectiveness of the fuzzy controller for controlling the drive system 
is demonstrated by comparing system response using the fuzzy controller, with that using a 
conventional PI controller. 

The computer program for simulation of proposed scheme is done in C language. Tlie 
Runge-Kutta fourth order method has been used for solving the differential equations 
governing system performance. 

2.2 Induction machine model 

2.2.1 Macliine voltage equations in synchronously rotating reference frame variables 

Equation of induction machine in two phase variables [8J can be written as 

Vqdos — TsiqdOs We^dqs"!" P^dOs (2.1) 

VqdOr~ rriqdOr’t' (We"Wj-)Xdqr +pA.qdOr (2-2) 

where 

(Xdq/=[Xds -Xq, OJ 

= (Xdr -Xqr 0) 

(VqdO.0‘ = [Vqs Vds 0] 

(VqdOr)^ =a (Vqr Vdr 01 

(iqdOs) “ fi<|S ids 0] 
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• X • 

(iqdOr) “ [iqr idr 0] 
and relation between flux linkages and currents 


^qdOs 


KsLs(Ks)-' 

KsLsr(K,)-' 


iqdOs 

^qdOr 


Kr(Lsr)''(Ks)-' 

KrUKr)-^ 


iqdOr 


The matrix KsLs(Ks)'’ can be expressed as 


KsLs(Ks)- 


and 


Lis+M 0 0 

0 Lis+M 0 

0 0 Lis 


Ks(Lsr)''(Ks)-’=KsLsr(Kr)-' = 


MOO 
0 M 0 
0 0 0 


where Ks is given by 


Ks = 


cos 0e 
sin 0e 

1 

2 


cos (0e-12O°) cos (0e+12O°) 
sin(0e-12O®) sin (0e+12O°) 


1 

2 


2 


After simplifying eq. 2.1 and eq. 2.2 the voltage equations can be written as 

Vqs ~ Ts^qs We X.(js P^s 


(2.3) 
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Vds = rids - WeA^,s + pXds (2.4) 

Vqr = TriqrH- (We-Wr)Xdr+ pXqr (2.5) 

Vdr = fridr- (We-Wr)Xqr+ pXdr (2.6) 

The equations for ilux linkages in expanded form 

As|s = (Lls+ M)iqs + Miqr (2.7) 

7-ds = (Lis+ M)ids + Midr . (2.8) 

Aqr = (Llr+ M)iqr + Miqs (2.9) 

Adr = (L]r+ M)idr + Mids (2. 10) 


2.2.2 Torque equation in two-phase variables 

The torque equation in two-phase variables can be written as 

Te=(P/2)[(Ks)-‘iqd0.j'"^ [Lsr](Ks)-'iqdOr (2.11) 

ao 

The above equation in terms of two phase stator flux and stator currents 
Te = (3/2)(P/2)fXdJq«-Xqdds) (2. 12) 


2.2.3 Calculation for slip speed 

For controlling the rotor flux of induction motor it is desirable that 


Xqr=0 (2.13) 

Xt = ./”(Adr^ + ) = Constant (2. 14) 

Since rotor is squirrel cage 

Vdr=Vqr=0 (2.15) 

From equations (2.5), (2.6), (2.9), (2. 10), (2. 15) 

idr=0 (2.16) 

(Wc-Wf). = slip speed = Wd = (-rAdr)iqr (2. 17) 

Xdr=Mids (2.18) 

V=-(M/Lr)iqs (2.19) 

combining equations (2.17,2.18,2.19) 

Wd = (r/LO(iqyids) (2.20) 
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where iq., is the torque component of cun'ent and ia., is the flux component of cuiTent. 

2.2.4 Calculation ol'Tiiree phase Reference Currents 

Calculation of three phase reference currents (ias*, ibs*, ics* ) involves the transformation of 
two pha.se reference currents I'rom synchronously rotating reference fiame to three phase 
stationary reference frame using a vector rotator (Fig 2. 1). 

The equation governing the transformation from two phase to tliree phase 


labcs 


Ka 


iqdOs 


( 2 . 21 ) 



from the Fig. 2. 1 

ias* = iqs* cos (9O+0e) + i<Js*COS (Ge) 
or, 

ias*= imCOS (0e+Y) (2-22) 

Similarly, 

ibs* = im COS (0e+ Y + 120°) (2.23) 

ics* = i„,cos(0e+Y-12O°) (2.24) 

where ' 

im=/"[(vy+(idsy] 

Y = tan'‘(iq5*/ids‘) 

The reference currents are fed to the bang-bang controller for controlling inverter switches. 
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2.3 Control sclienie 

Functional block diagram of indirect-field oriented control is given in Fig. 2.3. The control 
scheme basically concentrates on switching the PWM inverter so as to keep rotor flux of 
induction motor constant. The constant flux operation ensures independent control over 
torque of induction motor similar to torque control obtained from u separately excited fully- 
compensated dc motor. As in case of dc motor there are two orthogonal mmfs, one 
contributes for torque production while other contributes for establishment of working flux. 
In case of the induction motor tliis condition is not apparent, so from the point of view of our 
understanding the stator current vector is divided mto two components; one is called flux 
component of current while other is called torque component of cument. Since rotor flux can 
not be changed with a much faster rate due to large rotor time constant, the flux component 
of cun'ent is kept constant while torque component of current is incremented or decremented 
based on the torque requirement. The closed-loop control requires automatic adjustment of 
torque so that drive has good speed response under transient as well as steady-state operation. 
In the present work the fuzzy controller sets the torque component of current reference 
depending on magnitude of speed en or and change in speed eiror. The inverter is switched in 
such a manner that nK)tor cuiTents follow the set cuirent references within a hysteresis band 
(bang-bang control). The bang-bang logic can be summarized as follows. 
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Block diagram of indirect field-oriented control scheme 
Fig. 2.3 









If ia, ib, ic are motor currents and ias*, ibs* and ics* are the current reference and h is the 
hysteresis band 


If ia < (ias* - h) 
If ia > (ias* + h) 
If ib < (ibs* - h) 
If ib > (ibs* + h) 
If ic < (ics* - h) 
If ic > (ics* + h) 


Sa = 1 and Sa* = 0 
Sa = 0 and Sa* = 1 
Sb = 1 and Sb* = 0 
Sb = 0 and Sb* = 1 
Sc = 1 and Sc* - 0 
Sc = 0 and Sc* = 1 


2.4 Fuzzy logic controller simulation 

Fuzzy logic control is utilized for plants with complex dynamics that can not be known 
precisely. The fuzzy logic control approach for induction motor drive system is very useful 
since exact mathematical model of the system is not required. The function of fuzzy logic 
controller is to transform linguistic control rules into control strategy based on expert 
knowledge. 

The fuzzy logic control can be divided into four main functional blocks namely [9] 

1 . Knowledge base 

2. Fuzzification 

3. Inference mechanism 

4. Defuzzification 



error change torque current increment 

Fig. 2.4 Functional block diagram of fuzzy logic control 
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The knowledge base is composed of a data-base and a rule base. The data base, consisting of 
input and output membership functions, provides information for the appropriate 
fuzzification operations, the inference mechanism and defuzzification. The rule base is made 
up of a set of linguistic rules relating the fuzzy input variables to the desired control actions. 
Fuzzification converts a crisp input signal, the error (e), and error change (Ae), into fuzzified 
signals that can be identified by the levels of membership in the fuzzy sets. The inference 
mechanism uses the collection of linguistic control rules to convert the input conditions into 
fuzzified output. 

Finally, defuzification converts the fuzzy outputs into crisp control signals, which in the 
proposed system are the changes in torque component of current (Aiqs*) that drive the 
induction motor. Design procedure for fuzzy logic controller is as follows. 

The fuzzification maps the error, e(k) = (Wrs - Wr), and the error change Ae(k) = e(k)-e(k-l), 
to linguistic labels of the fuzzy sets. The proposed system uses following linguistic labels. 

For speed error fe) 

NL Negative Large 

NM Negative Medium 

NS Negative Small 

Z Zero 

PS Positive Small 

PM Positive Medium 

PL Positive Large 


Membership function 



Fig. 2.5 Membership function for speed error 
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For speed error change Ae 

NL Negative Large 
NS Negative Small 
Z Zero 

PS Positive Small 
PL Positive Large 

Membership function 



Ae-^ 

Fig. 2.6 Membership function for change in speed error 
For output current change A 
NL Negative Large 
NS Negative Small 
Z Zero 
PS Positive Small 
PL Positive Large 

Membership function 



Fig. 2.7 Membership function for change in torque component of current 
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2.5 Control rules 

Ejiowledge base involves defining the rules represented as statements governing the 
relationship between input and ou^ut variables in terms of membership functions. The 
control rules are represented as a set of if then rules, for example if speed error is PL and 
error change is Negative Large then change in output is Positive Large. Table2.1 shows the 
various fuzzy control mles. 

The control rules are evaluated by an inference mechanism. During the rule evaluation step, 
the combination of selected mles in knowledge base is evaluated. There are several 
mechanism for implementation of inference mechanism, such as min-max algorithm, the 
correlation algorithm etc [9]. The present work uses Mamdani algorithm for drawing 
inference. 

The output generated by fuzzy controller, which is used to control the plant, must be crisp 
value according to fuzzy output firom inference mechanism. The defuzzification block does 
this job i.e. it transforms crisp output in accordance with fuzzy output generated by inference 
mechanism. Many defuzzification strategies are available [9], such as, the weight average 
criterion, the mean-max (mean of maximum) membership, and center-of-area method 
(Centroid method). The center of area method is used here for defuzification, which can be 
summarized as follows. 

Suppose speed error is -5 and change in speed error is -10, then e will be NS or Z and 
change in error will be NS. So the rules relating the error and error rate are 


If e is NS and Ae is NS then change in output current is NS (Fig. 2.7). 
If e is Z and Ae is NS then change in output current is PS (Fig. 2.8). 
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Table 2. 1 Control Rules 


If speed error is e 

and error change is A e 

Then increment in output 
current is A iqs* 

PL 


PL 

PM 

NL 

PL 

PM 

NS 

PL 

PM 

Z 

PS 

PM 

PS 

PS 

PM 

PL 

NS 

PS 

NL 

PL 

PS 

NS 

PS 

PS 

Z 

PS 

PS 

PS 

PS 

PS 

PL 

NS 

Z 

NL 

PS 

Z 

NS 

PS 

Z 

Z 

z 

Z 

PS 

NS 

Z 

PL 

NS 

NS 

NL 

PS 

NS 

NS 

NS 

NS 

Z 

NS 

NS 

PS 

NS 

NS 

PL 

NL 

NM 

NL 

PS 

NM 

NS 

NS 

NM 

Z 

NS 

NM 

PS 

NL 

NM 

PL 

NL 

NL 


NL 
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Z NS PS 



Fig.2.8b 

Center of ai'ea gives the defuzzified value of current increment as 
C12 = (CiAi + C2A2 )/(Ai + A 2 ) 


2.6 Computer simulation of the drive system 

This section gives flowchart of the computer program developed for simulating the 
performance of the induction motor drive. 
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The proposed system has been simulated in 266 MHz pentium pc in Linux environment. The 
simulation program is written in C and the flowchart is given in Fig. 2.9. The complete C- 
program listing is included in Appendix B. 

2.7 Simulation results 

The drive has been given a step change in speed command and a step change in torque 
command to demonstrate feasibility of proposed scheme under varying operating conditioiis 
of operation. Computer traces of rotor flux response under change in speed command and 
torque command are shown in Fig. 2.10. The effectiveness of fuzzy controller is validated by 
comparing the results with that obtained using a conventional PI controller. 

2.7.1 Rotor flux response under change in speed command 

Fig. 2.10a, 2.10c show rotor flux response and rotor speed response when a step change of 
-1000 rpm (speed reversal)is applied at t =1.0s. Initially induction motor is allowed to nan at 
zero speed by setting speed command to zero. At t = 0.1s speed command is set at +500 rpm 
and further at t = 1.0s speed command is set at -SOOrpm. From Fig. 2.10a, rotor flux is 
settling at its steady state value in 0.2s from its value at standstill again in the region of speed 
reversal a slight overshoot is seen in rotor flux, in tliis region macliine is operating in 
regenerative braking mode. Rotor speed response (shown in Fig. 2. 10c) shows that macliine 
settles at +500 rpm within 0.1s and from +500 rpm to -500rpm in 0.2s. Fig. 2.10b shows 
locus of rotor flux in d-q plane under change in speed command, the circular nature of rotor 
flux vector in d-q plane ensures the vector control operation of induction motor drive system 
in transient as well as in steady state. 

2.7.2 Rotor flux response under change in torque command 

Fig. 2. lOd, 2. lOe, 2. lOf show rotor flux response under change in torque command applied to 
the proposed system. Initially machine is allowed to run at +500 rpm by setting speed 
command to +500 rpm and torque command to zero. At t =1.0 s. a step of +5.0 N.iu is 
applied. From Fig. 2. lOd, it can be seen that magnitude of the rotor flux is constant in spite of 
applied torque. Tlie circular nature of rotor flux vector (shown in Fig. 2. lOe) in d-q plane 
ensures vector control operation of induction motor drive system Torque response (shown in 
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Fig. 2. 1 Of) shows that electrical torque of inacliuie settles at set command witliiii a few 
supply-cycle thereby showing excellent torque response of the proposed system. 

2.7.3 Comparison 

Torque response and speed response of the proposed system is compared with the results 
obtained from operation of proposed system with PI controller. Simulated results are taken 
under change in speed command from -t-500 rpm to -500 rpm. Fig. 2. 1 1 shows flux, torque 
and speed response when induction motor drive is controlled by a fuzzy controller, and Fig. 
2. 12 shows flux, torque and speed response when induction motor drive is controlled by a PI 
controller. As can be seen from rotor speed response (Fig. 2.11c and Fig. 2.12c) that under 
drive operation with fuzzy controller rotor settles from +500 rpm to -500 rpm within 0.2 s 
whereas this figure is approximately 1.0 s for PI controller. Further from torque response 
curve (Fig. 2.11b and Fig. 2.12b) it can be concluded that torque response of the fuzzy 
controlled drive system is very fast as compared to that of the conventional PI controller. 
This discussion demonstrates the effectiveness of fuzzy controller for controlling non-linear 
systems like induction motor. 

2.8 Conclusion 

The proposed indirect field oriented control scheme is studied by computer simulation. The 
simulated results show that fuzzy logic control can give fast torque response under transient 
conditions. The comparison of simulated results with that obtained from conventional PI 
controller shows that fuzzy logic control is very fast. 
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Timc(secs)-> Timc(sccs)-> 

(c) Rotor speed (0 rvlcctromagnctic torque 

Wavcl{»frn :i j),c with step chaugc* ol'-KKK) fpm at t-!.Os VVavcfoim d,e.f with step change of -KKXlrpin at 


Roior flux waveform under chanpc in speed and torque command 

Fig. 2.10 






Rocor specd{rpm)>-> _ Rotor nux{ub)- 
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(a) Rotor (lux 



(b) Elcciromugnctic torque 




(c) Rotor speed (d) Machine phase current 

Fig. 2,1 1 Simulated results obtained from proposed system with controller as fuzzy 
(Step change of -1000 rpm at t = 1.5s) 
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Timc(scc)-> 


Tiine(scc)-> 


(c)Rolor speed 


(d)Machine phase current 


Fig. 2.12 Simulated results obtained from the proposed system with PI controller 
(Step change of -1000 rpm at t=1.5s) 
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Chapter 3 

Hardware and Control Software Implementation 

3.1 Introduction 

The proposed system has been practically implemented with the help of a personal computer 
having a liigh performance data acquisition card (PCL-208). The data acquisition card 
provides a channel for communication between software and hardware. There are some 
distinct advantages offered by PC-based implementation such as user friendliness, ease of 
programming and availability of large memory for storage. Additionally PCs have standard 
architecture so a software written for a paiticular PC cun be run on other PC manufactured by 
some other firm. This type of flexibility is not present with micro-controller systems, where a 
software written for a particular micro-controller system may not be executed on other 
micro-controller systems. 

The block diagram of the proposed system is shown in Fig. 3.1. Various blocks such 
as error generator, fuzzy controller, torque current limiter, vector rotator, slip frequency 
calculator and pulse width modulation logic has been implemented in real time by writing a 
control program in C which is given m the Appendix C. 

3.2 Specification of data acquisition card 

Tlie specification of data acquisition card is as follows [12J. 

1. Switch selectable 16 single-ended or 8 differential analog input channels. 

2. An industrial standard 12 Bit successive approximation converter (ADC674) to convert 
analog inputs. Tlie maximum A/D sampling rate is 60 KHz in DMA mode. 

3. Switch selectable versatile analog input ranges. 

Bipolar: -\-l- 0.5 V, +/- 1 V, +!- 2.5 V, -t-/- 5 V, +/- 10 V. 

Unipolar. + 1 V, + 2 V, -i- 5 V, + 10 V. 
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4. Provide three A/D trigger modes: software trigger, programmable pacer trigger and 
external trigger pulse trigger. 

5. A/D converted data can be transfen'ed by program control, interrupt handler routine or 
DMA transfer. 

6. An INTEL 8254 Programmable Timer/ Counter provides pacer output (trigger pulse) at the 
rate of 2.5 MHz to 71 minutes/ pulse to the A/D. The timer time base is switch selectable 10 
MHz or 1 MHz. One 16-bit counter channel is reserved for user configuration applications. 

7. Two 12 bit monolithic multiplying D/A output channels. Output range of 0 to +5 V can be 
created by using the onboard -5 V reference. This precision reference is derived from the 
A/D converter reference. External AC or DC reference can also be used to generate other 
D/A output ranges. 

8. TTL/ DTL compatible 16 digital input and 16 digital output channels. 

3.3 Setting of the data acquisition card 

The .setting of card for present application is as follows: 

1. Base address is set at 0x300 and corresponding I/O range is from 300 to 30F hexa- 
decimal. 

2. Analog input and analog output are set in bipolar mode and range of output voltage 
is±5V. 

3. Analog input channels are set in single ended mode. 

4. The clock frequency of 8254 is set at frequency of 10 MHz. 

For inputting actual speed, three phase actual currents and speed reference to the software 
five single ended analog input channels have been used. Tlie actual currents are compared 
with the reference currents and the PWM logic is generated. Tlie PWM signals to the upper 
transistors of the inverter are taken out by means of tliree digital output ports, wliich give 
+5V output corresponding to ON position of switch. Tlie firing pulse to lower switch of a leg 
is generated by complementing puKse of upper switch. In real time the complement of pulse 
of upper switch is obtained by using NOT gates. The use of NOT gates also ensures 
protection against any programming mistake, wliich may simultaneously ON upper and 
lower switch of a leg. 
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3.4 Lockout circuit 

Power devices require a finite amount of time to turn-off and tuni-on. At the time of tuni-oii 
of upper switch and tura-off of lower switch of a leg, there exist an interval in winch both 
tleviccs conduct. This type of situation should not exist in real time application since 
simultaneous conduction of both switches of a leg will short circuit the DC link voltage 
giving rise to high current thereby damaging power devices. To avoid tins situation a 
sufficient amount of blanking time should be provided between turn-on of upper switch and 
tuni-off of lower switch. In real time tins blanking time has been provided by using 74123 
itiono-stable multivibrator, which generates proper time delay governed by externally 
connected resistance and capacitance. In the proposed scheme a lock out time of 20 ji. s has 
i>cen chosen. The circuit diagram of lock out circuit is shown in Fig. 3.2. 

3-5 IGBT gate drive circuit 

Power devices require large amount of gate currents for turning them ON and OFF, tins 
requirement can not be met directly by outputs generated by low power control circuit, which 
provides pulses of 0-5V. The large gate current requirement is met by amplifying outputs of 
low power control circuitry. A pulse amplifier, wliich generates + 12V output corresponding 
to 0-5V generated by control circuit, has been used for experiment. Pulse amplifier can be 
designed to detect any over current present in the device by monitoring collector voltage. 
Plring circuit (shown in Fig. 3.3) used here gives over current indication by lighting up an 

led. 

Motor current sensing 

Actual current is sensed by Hall effect current sensor, which gives output at its measuring 
terminals in proportion with actual current. Use of Hall effect current sensor for current 
sensing also ensures isolation between high power circuit and low power control circuit. 
Output of cun'ent sensor used in experimental work, is .63V/amp. The detailed specification 
of the current sensor is given in Appendix A. 
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3.7 Control software 

The control software gives coniniand to the processor to read input data fironi specified 
hardware port and send certain output data to the specified hardware port. In the present 
a[iplicalion contiol .software rctids five input data I'roiu five analog input cliannels. These 
input data are used for determining the status of switching signals and finally these firing 
signals are outputted to specified hardware port. The execution time for reading five inputs 
ft-t)m hardware port and sending three outputs to hardware port is 90 micro seconds and clock 
speed of pentium processor is 266 MHz. The flowchart for deriving control program is given 
in Fig. 3.4. The control program is listed in Appendix C. 

3.8 Conclusion 

The proposed system is implemented on a tliree phase induction motor (specifications are 
given in Appendix A). The experimental results obtahied confirm the feasibility of the 
proposed system. 
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A view of lock out circuit for tlic leg of phase a 
Fig. 3.2 
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IGBT Gate Drive Circuit 


Start 



Fig. 3.4 Flow-chart for the control program 
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Chapter 4 

Experimental and Simulation Results 


4.1 Introduction 

The proposed scheme has been tested on a 3-phase, 415 V star connected induction motor, 
which is mechanically coupled with a 230 V, 1500 ipm dc generator for loading the 
induction motor indirectly. The magnitude for various parameters for inverter, motor and the 
dc generator are given in Appendix A. 

4.2 Speed response 

The proposed system is given a step change in speed ti'om +500 rpm to -500 tpm. Simulated 
speed response, rotor flux response and torque response during speed reversal is shown in 
Fig. 4.1. For t < 0.1 s machine is allowed to mn at zero speed by setting speed reference to 
zero, at t = 0. 1 s reference speed is set at +500 rpm, further at t = 1.0 s reference speed is set 
at -500 rpm. From simulated speed response it can be seen that machine settles at +500 rpm 
from standstill in approximately 0. 1 s, and at -500 rpm from +500 rpm witliiii 0.2 s. From t = 
1,0 s to t = 1.2 s machuie is operating in regenerative braking mode. From rotor flux response 
it can be seen that rotor flux settles in steady state value withiit 0.2 s, during speed reversal 
rotor flux has slight overshoot and after reversal rotor flux again settles at its steady state 
value. From the torque response curve, it can be seen that drive has very fast torque rise or 
fall during transient conditions, pulsations present in torque is inlierent in pulse width 
modulated inverter fed drives. Speed response and macliine current obtained from 
experiment are shown in Fig. 4.2 and Fig. 4.3, From which it can be seen that they are in 
clo.se agreement with the simulated speed and cur|ent responses. 

4.3 Torque response 

The induction motor in proposed system is subjected to a step change in torque command of 
+5 N.m at t = 1.0 s, from the simulated torque response (shown in Fig. 4.4) it can be seen that 
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torque resixinse is very fust. Rotor flux response shows that magnitude of rotor flux is 
constant under transient as well as in steady state conditions. Speed response shows that 
sudden application of torque causes a slight dip in speed but witliin a few time motor again 
settles at set speed. In the experiment the load test is performed on the induction motor by 
driving a dc generator feeding a resistive load, mechanically coupled to the induction motor. 
Tlie applied torque to the machuie is measured by measuruig armature current of dc 
generator. In the proposed scheme applied torque is 3.37 N.m. Experimental speed response 
shows that machine follows the set speed response after application of sudden load. 
Simulated speed response and practical speed response (shown in Fig. 4.5) ai'e found to be in 
close agreement. 

4.4 Current response 

Simulated current response and practical current response are shown in Fig. 4.6 and Fig. 4.7. 
From the simulated and practical cument responses it can be seen that actual currents are 
following the set current reference within a set hystcrisis band. The phase voltage waveform 
of the induction motor is shown in Fig. 4.8. 

4.5 Conclusion 

From the simulated and practical results obtained, following conclusions can be drawn. 

1. The induction motor is operating under constant rotor flux in steady state as well as in 
transient conditions, such as sudden application of step torque or application of step 
speed command. 

2. The constancy of rotor flux during transient and steady state indicates that macliine 
remains under vector control in transient as well as in steady state conditions. 
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(a)Ro[or speed during speed reversal from +500rpm to -500 rpm 



;i 0 0 r la'O i 


-- SOGrjsm 


Ti me csj i> 


(b) Ro(or speed during speed reversal from -i-500rpm to -500 rpm 

Sensitivity of. speed sensor: •04[)6v/rj)ni 

Rotor speed waveforms : (a) Simulated 

(b) Experimental 

Fig. 4.2 
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(a)Rt)tor speed 






(c) rotor speed and phase ‘a’ current waveforms 
scale :5.0v/cliv fur upper waveform 
scale: 2.0v/div for lower waveform 


Speed and curient wavefornis : a,b- Simulated 

c-l:xpcri mental 
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(c)Rlectromagnctic torque 
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Fig. 4.4 
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(a)phase ’a' reference and actual currents 



tfirDeCS) > 


(b)phasc ’a’ reference and actual currents (sacle:1.0v/div) 
sensitivity of current sensor 0.63v/atnp 

Reference and actual currents in phase ’a’ 

(a) Simulated 

(b) Experimental 

Pig. 4.6 
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currchi. 
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Fig. 4. 9 A view of practical setup 



Chapter 5 


Conclusion 


5.1 Salient features of the present work 
Salient features of present thesis are 

1 . Design and development of a fuzzy logic controller for the closed loop speed control of 
vector controlled induction motor drive. 

2. PC-based implementation of the vector controlled drive system. 

3. Comparison of simulation and experimental results. 

5.2 Scope for further research work 

The proposed scheme combines advantages offered by fiizzy logic controller for controlling 
induction motor drive system without precise knowledge of system behavior. But very 
precise control of torque and speed can not be achieved since control rules defining behavior 
of the plant, are itself fuzzy. Fuzzy logic controller can give very fast torque response in 
transient state as demonstrated by practical and simulation results. So as discussed in [10], 
good features of fuzzy logic controller in transient state and excellent speed regulation of 
phase locked loop can be combined to improve over all accuracy of the plant. 

In the proposed scheme the current references have been generated by calculating slip 
speed and adding it to actual rotor speed. The slip speed calculation is based on measunng 
rotor resistance and rotor inductance. Since resistance increases with increase in temperature, 
also with change in operating frequency the inductance changes. So with change in resistance 
and inductance, slip calculation no longer remains valid, this introduces error between actual 
slip speed and estimated slip speed. To account for variation in resistance due to temperature 
and frequency, calculation of resistance and inductance with varying operating conditions are 
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highly desiiable. So fuzzy control with rotor parameter adaptation is desirable under varying 
paiametric conditions. 

In the proposed scheme synchronous speed is calculated by sensing actual macliine speed in 
the form of tachometer output. RenK)vaI of speed seiivsor requires calculation of actual speed 
from known parameters such as voltage, cun'ents and constants of induction motor. A closed 
loop sensor- less fuzzy vector controlled drive can be an area of further investigation. 
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Spcdli cations of llie 3-pliasc, 50 Hz induciion motor 
Ratings: 

Power = 1.1 kW 
Voltage = 415 V line to line 
Cuirent = 2.6 A. 

Speed =1410 rpni 

9 

Moment of inertia of the motor-dc generator system = 0.01 kg m 
Rotor parameters: 

Resistance = 4.3 ohm 
Leakage inductance = 0.026 H 
Stator parameters: 

Resistance = 8.2 ohm 
Leakage inductance = 0.026 H 
Magnetizing inductance = 0.377 H 
Specifications of the Tachogenerator: 

Calibration factor = 6 V/ 1000 rpm. 

Resistance = 20 ohms. 

Max rpm = 2000. 

Specifications of dc-shunt generator coupled with the EM 
Ratings: 

Power ss 1 hp 
Voltage = 230 V. 
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CuiTent = 4.5 A. 

Speed = 1500 ipni 

I (i I Vr sp eci n cil t i « n.s : 

Voltage rating = 1200 V 
Cun-ent rating = 50 A 

Lock out delay time = 20 p s. 

Current sensor specifications 

The LEM module LA 55-P is a cuireiit transducer for electronic measurement of cuments 
with galvanic isolation between the primary (high power) and secondary (low power) 
circuits. 

sensitivity = .63 V/A 


Parameters for simulation 


DC link voltage 

350 V 

Hysterisis band gap 

±0.01 A 

Stator flux component of cuiTcnt 

2.0A 

Step size 

50 p s 

Sampling time for controller 

500 p s 

PI controller gains 

Kp = 5 


Ki = 50 

Parameters for hardware implementation 


DC link voltage 

350 V 

Hysterisis band gap 

±0.01 A 

Stator flux component of cuirent 

2.0 A 

Step size 

90 p s 

Sampling time for Controller 

900 p s 
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Listing of simulation program 


#include<stdio .h> 

# inc lude<ma th . h> 

#define PI 3.14159 

/* defining linguistic variables */ 


#def ine 

NLS 



e<=-150 

#def ine 

NMS 

e>=-200 

&5c 

e<=-40 

#def ine 

NSS 

e>=-50 

8 c&l 

e<=0 

#def ine 

ZS 

fD 

V 

ii 

1 

o 

ScSc 

e<=10 

#def ine 

PSS 

e>=0 


e<=50 

((define 

PMS 

e> = 40 

6rSc 

e<=200 

#def ine 

PLS 

e>=150 



#def ine 

NLC 



de<=-30 

#def ine 

NSC 

de>=-40 

ScSc 

de<=0 

#def ine 

ZC 

de>=-5 

ScSc 

de<=5 

#def ine 

PSC 

de>=0 

&& 

de<=40 

#def ine 

PLC 

de>=3 0 




/* Function for defuzzification */ 


float status (float e, float de) ; 

float z,u,v,x[30],er, ch, mship, a [30 ] ; 

float error ( float , int) ; 

float change (float, int) ; 

float min ( float , float ) ; 

int l,jl,kl; 

for (1=0;1<27;++1) 

{ 

a[l]=0; 

x[l]=0; 

} 

u=0 ; 
v=0 ; 

/* RULE 1 */ 

if(NLS) /*PL*/ 

{ 

jl=l; 
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niGhip=error ( e , j 1 ) ; 
a [0 ] = . 7 5*mship* (2-inship) ; 
x[0]=3 .5; 

} 

/* RULE 2, 3,4,5, 6 */ 

if (NMS) 

{ jl=2; 
if (NLC) 

{kl=l; 

er=error ( e , j 1 ) ; 
ch= change (de, kl) ; 
mship=min(er, ch) ; 
a [1] = . 75*mship* (2-mship) ; 
x[l]=3 .5; 

} 

if (NSC) 

{kl=2; 

er=error (e , j 1 ) ; 
ch=change (de, kl) ; 
inship=min (er , ch) ; 
a [2 ] = . 75*mship* (2-mship) ; 
x[2]=3.5; 

} 

if (ZC) 

{kl = 3; 

er=error (e, jl) ; 
ch=change (de, kl) ; 
mship=min (er , ch) ; 
a [3 ] =1 . 5*mship* (2-mship) ; 
'x[3]=1.5; 

} 

if (PSC) 

{kl=4; 

er=error (e, jl) ; 
ch=change (de, kl) ; 
mship=min (er , ch) ; 
a [4] =1. 5*mship* (2-mship) ; 
X [ 4 ] = 1 . 5 ; 

} 

if (PLC) 

{kl=5; 

er=error ( e, j 1) ; 
ch=change (de, kl) ; 
mship=min’(er , ch) ; 
a [ 5] =1 . 5*mship* (2-mship) ; 
x[5]=-1.5; 

} 


/*PL*/ 


/*PL*/ 


/*PS*/ 


/*PS*/ 


/*NS*/ 
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} 

/* RULE 7 8 9 10 11 */ 

if (NSS) 

{jl=3; 
if (NLC) 

{kl=l; 

er=error (e, jl) ; 
ch=change (de, kl) ; 
mship=min (er , ch) ; 
a [ 6 ] = . 75 *inship* (2 -mship) ; 
x[6]=3.5; 

} 

if (NSC) 

{kl=2; 

er=error (e, jl) ; 
ch=change (de , kl) ; 
mship=min ( er , ch) ; 
a [ 7 ] =1 . 5*mship* (2 -mship) ; 
x[7]=1.5; 

} 

if (ZC) 

{kl=3; 

er=error (e, jl) ; 
ch=change ( de , kl ) ; 
mship=min (er , ch) ; 
a [ 8] =1 . 5*mship* (2 -mship) ; 
x[8]=1.5; 

} 

if (PSC) 

{kl=4; 

er=error (e, jl) ; 
ch=change (de, kl) ; 
mship=min(er, ch) ; 
a [ 9 ] = 1 . 5 *mship * ( 2 -mship ) ; 
x[9]=1.5; 

} 

if (PLC) 

{kl=5; 

er=error (e, jl) 
ch=change (de , kl) ; 
mship=min ( er , ch) ; 
a [ 10 ] =1 . 5 *mship* ( 2 - mship) ; 
x[10j =-l. 5; 

} 

} 

/* RULE 12 13 14 15 16 */ 


/*PL*/ 


/*PS*/ 


/*PS*/ 


/*PS*/ 


/*NS*/ 
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if (ZS) 

{ jl = 4; 
if (NLC) 

{kl=l; 

er=error ( e , j 1 ) ; 
ch=change (de, kl) ; 
mship=min (er , ch) ; 
a [11] =1 . 5*mship* (2-mship) ; 
x[ll]=1.5; 

} 

if (NSC) 

{kl=2; 

er=error (e, jl); 
ch=change (de, kl) ; 
mship=min (er, ch) ; 
a [ 12 ] =1 . 5 *mship* (2-mship) ; 
x[12]=1.5; 

} 

if (ZC) 

{kl=3; 

er =error ( e , j 1 ) ; 
ch=change (de, kl) ; 
mship=min ( er , ch) ; 
a [13 ] = . 5*mship* (2-mship) ; 
x[13]=0; 

} 

if (PSC) 

{kl=4; 

er=error (e, jl); 
ch=change (de, kl) ; 
mship=min ( er , ch) ; 
a [14] =1 . 5*mship* (2-mship) ; 
x[14]=-1.5; 

} 

if (PLC) 

{kl=5; 

er=error (e, jl) ; 
ch=change ( de , kl ) ; 
mship=min (er , ch) ; 
a [ 15 ] =1 . 5*mship* (2-mship) ; 
x[15]=-1.5; 

} 

} 

/* RULE 17,18,19,20,21 */ 

if (PSS) 

{jl=5; 
if (NLC) 


/*PS*/ 


/*PS*/ 


1*7,* I 


/*NS*/ 


/*NS*/ 


/*PS*/ 
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{kl=l; 

er=error ( e , j 1 ) ; 
ch=change (de, kl) ; 
mship=min (er , ch) ; 
a [ 16 ] =1 . 5*mship* (2-mship) ; 
X [16] =1.5; 

) 

if (NSC) 

{kl=2; 

er=error (e, jl) ; 
ch=change (de, kl) ; 
mship=min (er , ch) ; 
a [ 17 ] =1 . 5*mship* (2-mship) ; 
x[17]=-1.5; 

} 

if (ZC) 

{kl=3; 

er=error ( e, j 1) ; 
ch=change (de, kl) ; 
mship=min (er, ch) ; 
a [ 18 ] =1 . 5 *mship* ( 2 -mship) ; 
x[18]=-1.5; 

} 

if (PSC) 

{kl=4; 

er=error (e, jl) ; 
ch=change (de, kl) ; 
mship=min (er , ch) ; 
a [ 19 ] =1 . 5*mship* (2-mship) ; 
x[19]=-1.5; 

} 

if (PLC) 

{kl=5; 

er=error ( e , j 1 ) ; 
ch=change (de , kl) ; 
mship=min (er , ch) ; 
a [20 ] = . 75*mship* (2-mship) ; 
x[20]=-3.5; 

} 

} 

/* RULE 22,23,24,25,26 */ 

i f ( PMS ) 

{ jl=6; 
if (NLC) 

{ klssl ; 

er=error (e, jl) ; 
ch=change ( de , kl ) ; 


/*NS*/ 


/*NS*/ 


/*NS*/ 


/*NL*/ 


/*PS*/ 
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inship=min ( er , ch) ; 

a [21] =1 . 5*mship* (2-mship) ; 

x[21]=1.5; 

} 

if (NSC) 

{kl=2; 

e r ~ c r r o r ( o , j 1 ) ; 
ch=change (de, kl) ; 
mship=min (er, ch) ; 
a [22 ] =1 . 5*mship* (2-mship) ; 
x[22]=-1.5; 

} 

if (ZC) 

{kl=3; 

er=error ( e , j 1 ) ; 
ch=change (de, kl) ; 
mship=min (er , ch) ; 
a [23 ] =1 . 5*mship* (2-mship) ; 
x[23]=-1.5; 

} 

if (PSC) 

{kl=:4; 

er=error (e, jl) ; 
ch=change ( de , kl ) ; 
mship=min (er , ch) ; 
a [24 ] = . 75*mship* (2-mship) ; 
x[24]=-3 .5; 

} 

if (PLC) 

{kl=5; 

-er=error (e, jl) ; 
ch=change (de , kl ) ; 
mship=min ( er , ch) ; 
a [25] =. 75*mship* (2-mship) ; 
x[25]=-3 .5; 

} 

} 

/* RULE 27 */ 

if (PLS) 

{jl=7; 

mship=error (e,jl); 

a [2 6 ] = . 7 5*mship* ( 2 -mship) ; 

x[26]=-3.5; 

} 

f or ( 1=0 ; 1<27 ; ++1) 

{ 

u+=a[l]*x[l] ; 


/*NS*/ 


/*NS*/ 


/*NL*/ 


/*NL*/ 


/* NL*/ 
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v+=a [1] ; 

a[l]=0; 

x[l]=0; 

} 

z=u/v; 
return ( z ) ; 

} 

/* FUNCTION FOR SPEED ERROR */ 

float error (float e,int jl) 

{ 

float g; 

/* NEGATIVE LARGE */ 

switch ( j 1 ) 

{ 

case 1: 
if (NLS) 

{ 

if (e>=-300)g=-(l+e/150) ; 
else g=l; 

) 

break; 

/* NEGATIVE MEDIUM */ 
case 2 : 
if (NMS) 

{ 

if (e<=-120)g=(e+200) /80; 
else g=- (e+40) /80; 

} 

break; 

/* NEGATIVE SMALL */ 

case 3 : 
if (NSS) 

{ 

if (e<=-25)g=(e+50) /25; 
else g=-e/25; 

} 

break; 

/* ZERO */ 

case 4 : 
if (ZS) 

{ 

if (e<=0) g= (e+10) /lO ; 
else 

g=(10-e) /lO; 

} 

break; 
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/* POSITIVE SMALL */ 
case 5: 
if (PSS) 

{ 

if (e<=25)g=e/25; 
else g=(50-e)/25; 

) 

break; 

/* POSITIVE MEDIUM */ 
case 6 : 
if (PMS) 

{ 

if (e<=120)g=(e-40) /SO; 
else g= (200-e) /SO ; 

} 

break; 

/* POSITIVE LARGE */ 
case 7 : 
if (PLS) 

{ 

if (e<=300) g=e/150-l; 

) 

break ; 

} 

return (g) ; 

} 

/* FUNCTION FOR ERROR CHANGE */ 

float change (float de,int kl) 

{ 

float k; 
switch (kl) 

{ 

/* NEGATIVE LARGE */ 

case 1: 
if (NLC) 

{ 

if (de>=-60)k=-(de+30) /30; 
else k=l; 

} 

break; 

/* NEGATIVE SMALL */ 

case 2 : 
if (NSC) 

{ 
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if {de<=-20) k= (de+40) /20; 
else k=-de/20;; 

} 

break; 

/♦ZERO*/ 
case 3 : 
i f (ZC) 

{ 

if (de< = 0) k= (de+5) /5; 
else 

k= (5-de) /5; 

} 

break; 

/* POSITIVE SMALL */ 
case 4 : 
if (PSC) 

{ 

if (de<=20) k=de/20; 
else 

k=(40-de) /20; 

} 

break; 

/* POSTIVE LARGE */ 
case 5 : 
if (PLC) 

{ 

if (de<=60)k=(de-30) /30; 
else k=l; 

} 

break; 

} 

return (k) ; 

} 

/* FUNCTION FOR MINIMUM MEMBERSHIP */ 

float inin(float me, float mde) 

{ 

if (me>mde) return (me) ; 
else return (mde); 

} 

main ( ) 

{ 

int count=-l , lo=0 , list , po=4 ; 
float 

el=0 , eO , wrs , iqss=0 , idss , Kp, Ki , THsl=0 , iqsm, THs2 , THr=0 , THe, h 
float sir, Vdc, Im, ias, ibs, ics; 

float iqso, e, de, t0=0 , T,m[5] , n [5 ] , o [5 ] , p [5] , g, ia, ib, ic; 
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float M= .357, Llr= .026,Lls=.026,d,f=60. ,t=0, va , vb, vc ; 

float wr so , vds , vqs , iqs , iqr , ids , idr , J= . 0 12 ; 

float 

rs=8.25,r=4.32,_t, siqr=0 , siqs=0 , sidr=0 , sids=0 , Te, Tl=0 , wr=0 , B=( 


float 

s iqr 1=0 , s iqs 1=0 , s idr 1=0 , s ids 1=0 , iqsl=0 , iqr 1=0 , ids 1=0 , idr 1=0 ; 
float vn, vao, vbo, vco, z; 
char ch; 

FILE *fp ,*fn ,*fm; 

Fp = fopen("r" ,"w"); 

Fm = fopen("t", "w"); 

Fn = fopen("s", "r"); 

printf (" Enter the controller you want to choose Fuzzy or PI P 
,F\n") ; 

ch=toupper (getchar ( ) ) ; 

print f (" Enter the value h vdc kp ki idss iqsm T wrs") ; 

fscanf( fn, "%f %f %f %f %f %f %f %f %f \n" , &h, &_t , &Vdc, &Kp, 

ScKi, Scidss, Sciqsm, ScT, &wrs); 

g=M* (Llr+Lls) +Llr*Lls; 

wrso=wrs ; 

do 

{++count; 

if ( t< . 01) wrs=0 ; 

else if(t>=.01)wrs=wrso; 

if ( t> . 3 ) wrs=-104 . 72 ; 

if (count==10) {count=0; 

el=-wrs+wr ; 

e=el ; 

de= (el-eO) ; 
switch ( ch) 

{ 

case 'P' : 

iqss+=Kp* (-el+eO) -Ki*e0*T; 
break ; 
case ' F ' : 

iqss+= . 07*status (10*e, 1000*de) ; 
break; 

} 

e 0 = e 1 ; 

if (iqss>iqsm) iqss=iqsm; 
if (iqss<-iqsm) iqss=-iqsm; 

THsl+= (iqss/idss) * (r/ (Llr+M) ) *T; 

THs2=atan2 ( iqss , idss) ; 

THr+=wr*T; 

THe=THsl+THs2+THr; 

Iin=sqrt (iqss*iqss+idss*idss) ; 
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/* calculating current references */ 

ias=Im*cos (THe) ; 
ibs=Im*cos (THe-2*PI/3) ; 
ics=Im*cos (THe+2*PI/3) ; 

} 

/* bang-bang logic */ 

if ( ia< { ias-h) ) vao=Vdc ; 
if ( ia> ( ias+h) ) vao=0 ; 
if ( ib< ( ibs-h) ) vbo=Vdc ; 
if ( ib> ( ibs+h) ) vbo=0 ; 
if ( ic< ( ics-h) ) vco=Vdc ; 
if {ic> (ics+h) ) vco=0; 
vn= (vao+vbo+vco) /3 . ; 
va=vao-vn; 
vb=vbo-vn; 
vc=vco-vn; 

vds= (va- . 5* (vb+vc) ) *2 . /3 ; 
vqs= ( . 866* (vb-vc) ) *2 . /3 ; 

/* Runge-Kutta fourth order method */ 
f or ( list=0 ; list<4 ; ++list ) 

{ 

iqs= ( (Llr+M) *siqs-M*siqr) /g,' 
iqr= ( (Lls+M) *siqr-M*siqs) /g; 
ids= ( (Llr+M) *sids-M*sidr) /g; 
idr= ( (Lls+M) *sidr-M*sids) /g; 

m[list] = (vqs-rs*iqs) *_t; 
n [ list ] = (vds-rs*ids) *_t ; 
o [list] = (-r*iqr+wr*sidr) *_t; 
p [ list] = ( -r*idr-wr*siqr ) *_t; 
if ( list ! =3 ) 

{ 

siqs+=m[list] ; 
sids+=n[list] ; 
siqr+=o [list] ; 
sidr+=p [list] ; 

} 

) 

iqsl= ( (Llr+M) *siqsl-M*siqrl) /g; 
iqrl= ( (Lls+M) *siqrl-M*siqsl) /g; 
idsl= ( (Llr+M) *sidsl-M*sidrl) /g; 
idrl= ( (Lls+M) *sidrl-M*sidsl) /g; 
siqsl+= (m[0 ] +2 *m[l] +2*m[2 ] +in[3 ] ) / 6 . ; 
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sidsl+= (n[0] +2*n[l] +2*n[2] +n[3] ) / 6. ; 
siqrl+= (o [0] +2*o[l] +2*o[2] +o [3])/6.; 
sidrl+={p[0] +2*p[l]+2*p[2] +p[3] ) /6. ; 

Te= (3/2.)*(po/2)* (sidsl*iqsl-siqsl*idsl) ; 
wr+= (po/2 ) * ( 1/ J) * (Te-Tl-B* (2 . /po) *wr) *_t ; 
ia=idsl ; 

ib=" . 5* idsl+ . 866* iqsl ; 
ic=- . 5*idsl- . 866*iqsl ; 
t+=_t ; 

sir=sqrt ( sidrl*sidrl+siqrl*siqrl) ; 

/ * Output * / 

fprintf ( fp, " % . 5f %.4f %.4f %-4f %.4f %.4f %.4f %.4f 

% . 4f \n" , t , Tl , sidr , siqr , sir, ( 60/ (po*PI) ) *wr , Te, ia, ias+h) ; 

siqs=siqsl ; 

sids=sidsl ; 

siqr=siqrl ; 

sidr=sidrl ; 

} while ( t< . 6 ) ; 

} 
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Appendix C 

Listing of control program 


#include<stdio .h> 
include<conio . h> 
#include<dos . h> 
#include<time . h> 
#include<math. h> 

#define PI 4 . 0*atan ( 1 . 0 ) 


/* defining linguistic labels */ 


#def ine 

NLS 



e<=-150 

#def ine 

NMS 

e>=-200 

&& 

O 

1 

II 

V 

0) 

#def ine 

NSS 

e>=-50 

&& 

e<=0 

#def ine 

ZS 

(D 

V 

il 

1 

o 

&& 

(D 

A 

ii 

O 

#def ine 

PSS 

o 

ii 

A 

0) 

&& 

O 

in 

ii 

V 

0) 

#def ine 

PMS 

o 

II 

A 

0) 

ScSc 

e'^ — 200 

)fdef ine 

PLS 

e>=150 



#def ine 

NLC 



de<=-30 

#def ine 

NSC 

de>=-40 

ScSc 

de<=0 

ttdef ine 

ZC 

de>=-5 

ScSc 

de<=5 

#def ine 

PSC 

de>=0 

&& 

de<=40 

#def ine 

PLC 

de>=3 0 




double ainO , ainl , ain2 , ain3 , ain4 ; 
int port=0x300; 

/* function for defuzzification */ 

float status (float e, float de) 

{ 

float z,v,a[30] ,er, ch, mship, x[30 ] ; 

float error ( float , int) , u=0 ; 

float change (float, int) ; 

float min ( float , float) ; 

int l,jl,kl; 

for(l=0;l<27;++l) 

{ 

a[l]=0; 

x[l]=0; 

} 

u=0 ; 
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v=0 ; 

/* RULE 1 */ 

if (NLS) 

{ 

jl=l; 

inship=error ( e , j 1 ) ; 
a [ 0 ] = . 75 *mship* ( 2 mship) ; 
x[0] =3 . 5; 

} 

/ * RULE 2 , 3 , 4 , 5 , 6 * / 

if (NMS) 

{jl=2; 
if (NLC) 

{kl=l; 

er=error ( e , j 1 ) ; 
ch=change (de, kl) ; 
mship=min(er, ch) ; 
a [ 1 ] = . 75 *mship* ( 2 -mship) ; 
x[l]=3 .5; 

} 

if (NSC) 

{kl=2; 

er=error (e, jl); 
ch=change (de , kl ) ; 
mship=min(er, ch) ; 
a [ 2 ] = . 7 5 *mship* ( 2 -mship ) ; 
x[2]=3 .5; 

} 

if (ZC) 

{kl=3; 

er =error ( e , j 1 ) ; 
ch=change (de, kl) ; 
mship=min(er , ch) ; 
a [3 ] =1 . 5*mship* (2 -mship) ; 
x[3]=l. 5; 

} 

if (PSC) 

{kl=4; 

er=error ( e , j 1) ; 
ch=change ( de , kl ) ; 
mship=min(er, ch) ; 
a [4] =1 . 5*mship* (2-mship) ; 
x[4]=1.5; 

} 

if ( PLC ) 

{kl=5; 

er=error (e, jl) ; 


/*PL*/ 


/*PL*/ 


/*PL*/ 


/*PS*/ 


/*PS*/ 


/*NS*/ 
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ch=change (de, kl) ; 
mship=min (er , ch) ; 
a [5] =1 . 5*mship* (2-mship) ; 
x[5] =-l. 5; 

} 

} 

/* RULES 7,8,9,10, 11 */ 

if (NSS) 

{ D1=3; 
if (NLC) 

{kl=l; 

er=error (e, jl) ; 
ch=change (de, kl) ; 
mship=min (er , ch) ; 
a [ 6 ] = . 75*mship* (2-mship) ; 
x[6]=3.5; 

} 

if (NSC) 

{kl=2; 

er=error (e, j 1) ; 
ch=change ( de , kl ) ; 
mship=min (er , ch) ; 
a [7 ] =1 . 5*mship* (2-mship) ; 
x[7]=1.5; 

) 

if (ZC) 

{kl=3; 

er=error (e, jl); 
ch=change (de, kl) ; 
mship=min (er , ch) ; 
a [8] =1 . 5*mship* (2-mship) ; 
x[8I=1.5; 

} 

if (PSC) 

{kl=4; 

er=error (e, jl) ; 
ch=change (de, kl) ; 
mship=min (er , ch) ; 
a [ 9 ] =1 . 5*mship* (2-mship) ; 
x[9]=1.5; 

) 

if (PLC) 

{kl=5; 

er=error (e, jl); 
ch=change (de, kl) ; 
mship=min (er, ch) ; 
a[10}=:l. 5 *mship* (2-mship) ; 


/*PL*/ 


/*PS*/ 


/*PS*/ 


/*PS*/ 


/*NS*/ 
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x[10]=-1.5; 

} 

} 

/* RULES 12, 13. 14, 15, 16 

if (ZS) 

{ j 1=4 ; 
if (NLC) 

{kl=l; 

er=error (e, j 1) ; 
ch=change (de, kl) ; 
mship=min { er , ch) ; 
a [11] =1 . 5*mship* (2-mship) ; 
x[ll]=1.5; 

} 

if (NSC) 

{kl=2; 

er=error ( e , j 1 ) ; 
ch=change (de, kl) ; 
mship=min(er , ch) ; 
a [12 ] =1 . 5*mship* (2-mship) ; 
x[12]=1.5; 

} 

if (ZC) 

{kl=3; 

er=error (e, jl) ; 
ch=change (de, kl) ; 
mship=min (er , ch) ; 
a [13 ] = . 5*mship* (2-mship) ; 
x[13]=0; 

} 

if (PSC) 

{kl=4; 

er=error ( e , j 1 ) ; 
ch=change ( de , kl ) ; 
mship=min (er , ch) ; 
a [ 14 ] =1 . 5*mship* (2-mship) ; 
x[14]=-1.5; 

} 

if (PLC) 

{kl=5; 

er=error ( e , j 1 ) ; 
ch=change (de, kl) ; 
mship=min(er, ch) ; . 
a [15] =1 . 5*mship* (2-mship) ; 
x[15]=-1.5; 

} 


*/ 

/*PS*/ 

/*PS*/ 

I*Z*l 

/*NS*/ 

/*NS*/ 
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} 

/* RULES 17, 18, 19, 20, 21 */ 

if (PSS) 

{ jl=5; 

if(NLC) /*PS*/ 

{kl=l; 

er=error (e, j 1) ; 
ch=change (de, kl) ; 
mship=min ( er , ch) ; 
a[16] =1.5*mship* (2-mship) ; 
x[16]=1.5; 

} 

if (NSC) /*NS*/ 

{kl=2; 

er=error (e, j 1) ; 
ch=change ( de , kl ) ; 
mship=inin (er , ch) ; 
a [ 17 ] =1 . 5*mship* (2-mship) ; 
x[17]=-1.5; 

} 

if(ZC) /*NS*/ 

{kl=3; 

er=error (e, jl) ; 
ch=change (de, kl) ; 
mship=min(er, ch) ; 
a [ 18 ] =1 . 5*mship* (2-mship) ; 
x[18]=-1.5; 

} 

if(PSC) /*NS*/ 

{kl=4; 

er=error (e, j 1) ; 
ch=change (de, kl) ; 
mship=min (er , ch) ; 
a [19 ] =1 - 5*mship* (2-mship) ; 
x[19]=-1.5; 

} 

if(PLC) /*NL*/ 

{kl=5; 

er=error (e , j 1) ; 

ch=change (de, kl) ; 

mship=min (er , ch) ; 

a[2 0] = .75*mship* (2-mship) ,- 

x[20]=-3.5; 

} 

) 

/* RULES 22, 23, 24, 25, 26, */ 

if (PMS) 
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{jl=6; 
if (NLC) 

{kl=l; 

er =erroV ( e , j 1 ) ; 
ch=change (de, kl) ; 
mship=min (er , ch) ; 
a [ 2 1 ] = 1 . 5 *mship * ( 2 -mship ) ; 
x[21]=1.5; 

} 

if (NSC) 

{kl=2; 

er=error (e, jl); 
ch=change (de, kl) ; 
n\ship=min (er , ch) ; 
a [ 22 ] =1 . 5 *mship* (2 -mship) ; 
x(22]=-1.5; 

} 

if (ZC) 

{kl=3; 

er=error ( e , j 1 ) ; 
ch=change (de , kl) ; 
mship=min (er , ch) ; 
a [23 ] =1 . 5*mship* (2 -mship) ; 
x[23] =-l. 5; 

} 

if (PSC) 

{kl=4; 

er=error (e, jl) ; 
ch=change (de, kl) ; 
mship=min (er , ch) ; 
a[24] =. 75*mship* (2-mship) ; 
x[24] =-3 . 5 ; 

) 

if (PLC) 

{kl=5; 

er=error ( e , j 1 ) ; 

ch=change (de, kl) ; 

mship=min (er , ch) ; 

a[25] =. 75*mship* (2-mship) ; 

x[25]=-3.5; 

} 

) 

/* RULE 27 */ 

if (PLS) 

{jl=7; 

mship=error (e, j 1) ; 
a [ 2 6] = . 75*mship* (2-mship) ; 


/*PS*/ 


/*NS*/ 


/*NS*/ 


/*NL*/ 


/ *NL* / 


/* NL*/ 
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x[26]=-3.5; 

) 

for (1=0;1<27;++1) 

{ 

u+=a[l] *x[l] ; 
v+=a [1] ; 
a[l]=0; 
x[l]=0; 

} 

z=u/v; 
return ( z ) ; 

} 

/* FUNCTION FOR SPEED ERROR */ 

float error (float e,int jl) 

{ 

float g; 

/* NEGATIVE LARGE */ 

switch ( j 1) 

{ 

case 1: 
if (NLS) 

{ 

if (e>=-3 00)g=-(l+e/150) ; 
else g=l; 

} 

break; 

/* NEGATIVE MEDIUM */ 
case 2 : 
if (NMS) 

{ 

if (e<=-120)g=(e+200) /80; 
else g=- (e+40) /80; 

} 

break; 

/*, NEGATIVE SMALL */ 
case 3 : 
if (NSS) 

{ 

if (e<=-25) g= (e+50) /25; 
else g=-e/25; 

) 

break; 

/* ZERO */ 

case 4 : 
if (ZS) 

{ 
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if {e<=0)g=(e+10) /lO; 
else 

g=(10-e) /lO; 

} 

break; 

/* POSITIVE SMALL */ 
case 5: 
if (PSS) 

{ 

if {e<=25)g=e/25; 
else g= (50-e) /25; 

} 

break; 

/* POSITIVE MEDIUM */ 
case 6 : 
if (PMS) 

{ 

if (e<=120)g={e-40) /80; 
else g= (200-e) /80 ; 

} 

break; 

/* POSITIVE LARGE */ 
case 7 : 
if (PLS) 

{ 

if(e<=300)g=e/150-l; 

} 

break; 

} 

return (g) ; 

} 

/* FUNCTION FOR ERROR CHANGE */ 

float change (float de,int kl) 

{ 

float k; 
switch (kl) 

{ 

/* NEGATIVE LARGE */ 

case 1: 
if (NLC) 

{ 

if (de>=-60)k=- (de+30) /30; 
else k=l; 

} 
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break ; 

/* NEGATIVE SMALL */ 
case 2 : 
if (NSC) 

{ 

if {de<=-20) k= (de+40) /20; 
else k=- de/20 ; ; 

} 

break ; 

/* ZERO */ 
case 3 : 
if (ZC) 

{ 

if (de<=0) k= (de+5) /5; 
else 

k= (5-de) /5; 

} 

break ; 

/* POSITIVE SMALL */ 
case 4 : 
if (PSC) 

{ 

if (de<=20) k=de/20; 
else 

k=(40-de) /20; 

} 

break ; 

/* POSTIVE LARGE */ 
case 5 : 
if (PLC) 

{ 

if (de<=60)k=(de-30) /30; 
else k=l; 

} 

break; 

} 

return (k) ; 

} 

/* FUNCTION FOR MINIMUM MEMBERSHIP */ 

float ininCfloat me, float mde) 

{ 

if (me>mde) return (me) ; 
else return (mde) ; 

} 

main ( ) 

{ clock_t start, end; 
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int sa=0 , sb=0 , sc=0 ; 

int count=19; 

float el=0,e0=0,wrs=- 

104.72, iqss = 0 , idss=2 . 0 , Kp=5 . 0 , Ki=50 . 0 , THsl = 0 , iqsm=3 . 0 , THs2 = 0 , T 
Hr=0 , THe=0 , h=0 . 001; 

double lin=0 , ias = 0 , ibs = 0 , ics=0 , r=4 . 32 ; 

double e=0,de=0,T=. 001780, ia=0, ib=0 , ic=0 ; 

double M=.377,Llr=.026,d, t=0,- 

double _t= . 000089 ,wr=0 ; 

int chl=0 , ch2=0 , ch3=0 ; 

void fireout (int, int , int) ; 

void csensor ( void) ; 

char ch; 

FILE *fn; 

fn=fopen { "r " , "w" ) ; 

printf (" Enter the controller you want to choose Fuzzy or PI P 
,F\n" ) ; 

ch=toupper (getchar ( ) ) ; 

start=clock ( ) ; 
do 
{ 

++count ; 
csensor ( ) ; 

if (count==10 ) {count=0; 

wrs=32 *ain4 ; 

if ( t>l . 5 ) wrs=104 . 72 ; 

wr=34 . 89*ain2 ; 

el=-wrs+wr ; 

e=el ; 

de=el-e0 ; 

/* Fuzzy controller simulation */ 

switch (ch) 

{ 

case ' P ' ; 

iqss+= . 2 * (Kp* ( -el+eO ) -Ki*e0*T) ; 
case ' F ' : 

iqss+=0 . 07status (10*e, 1000*de) ; 
break; 

} 

e0=el ; 

if ( iqss>igsm) iqss=iqsm; 
if ( iqss<-iqsm) iqss=-iqsin; 

THsl+a (iqss/idss) * (r/ (Llr+M) ) *T; 

THs2=atan2 ( iqss , idss) ; 

THr+=wr*T; 
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THe=THs l+THs2 +THr ; 

Im=sqrt ( iqss*iqss+idss*idss) ; 

/* Setting current reference */ 
ias = Iin*cos (THe) ; 
ibs=Im*cos (THe-2*PI/3) ; 
ics = Im*cos (THe+2*PI/3) ; 

chl = floor(4095*(1.5-(1.8/2 8.6) *ias) /5.0) ; 

ch2 =chl/ 16 ; 

ch3 = 16* (chl-ch2*16) ; 

outport (port+4 , ch3 ) ; 

outport (port+5 , ch2 ) ; 

} 

ib=-l. 5713*ain0; 
ic=-l . 5713*ainl; 
ia=-l . 5713*ain3 ; 

/* bang-bang logic */ 
if ( ia< ( ias-h) ) sa=l ; 
if ( ia> { ias+h) ) sa=0 ; 
if ( ib< ( ibs-h) ) sb=l ; 
if ( ib> ( ibs+h) ) sb=0 ; 
if ( ic< { ics-h) ) sc=l; 
if ( ic> ( ics+h) ) sc=0 ; 
f ireout (sa, sb, sc) ; 
t+=_t ; 

} while (t<5 . 0) ; 
end=clock ( ) ; 

printf ( "T‘ime=%f T=%f " , (end-start) /CLK_TCK,T) ; 
outport (port+3 , 0x00 ); /* Inhibiting firing pulses */ 
} 

/* FUNCTION FOR SENDING FIRING PULSES TO OUT PORT 
void f ireout (int sa, int sb,int sc) 

{ 

int mp; 

mp=4*sc+2 *sb+sa; 
switch (mp) { 
case 0: 

outport (port+3 , 0x00 ) ; 

break; 

case 1: 

outport (port+3 , 0x01) ; 
break; 
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case 2 : 

outport (port+3 , 0x02 ) ; 
break ; 
case 3 : 

outport (port+3 , 0x03 ) ; 

break; 

case 4 : 

outport (port+3 , 0x04) ; 

break; 

case 5 : 

outport (port+3 , 0x05 ) ; 
break ; 
case 6 : 

outport (port+3 , 0x06) ; 

break; 

case 7 : 

outport (port+3 , 0x07 ) ; 
break ; 

) 

return; 

} 

/* FUNCTION FOR SENSING CURRENT */ 

void csensor (void) 

{ 

int st,uni=l, s_end=l, dif f ,bip, start=0, stop=4; 
int dtl, dth, adl, adt , c_reg, s_ch_val, r_ch; 
int chv,pool=0; 
char number ; 
int val; 


/* ******** step 1: INITIALIZE & SELECTE SOFTWARE TRIGGER 

★ *★*★★★* ifcr i»r ★ 

/* PORT+9 CONTROL REGISTER 

*/ 

val=0x70 ; 

outportb (port+9 , val) ; 
c_reg=inportb (port+9) ; 
if (c_reg != val) 

{ 

printf ("pcl-718 hardware verification failed! \n"); 
exit ( 0 ) ; 

} 

outportb (port+8,1); /* clear interrupt request */ 
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STEP 3: SET SCAN CHANNEL RANGE 


***************************** 


7 


S_ch_val=Stop*16+start; /* SET SCAN CHANNEL VALUE */ 

outportb (port+2 , s_ch_val) ; 

r_ch=inportb (port+2); /* READ BACK CHANNEL VALUE */ 

if (r_ch != s_ch_val) 

{ 

printf ("set scan channel f ailed ! \n" ) ; 
exit ( 0 ) ; 

} 

/* ******** step 4: PERFORM SINGLE A/D CONVERSION ********** 

*/ 

do { 

outportb (port,0); 
reread: st=inport (port+8 ) ; 

if ((st & 0x80)==0x80) 
goto reread; 

dtl=inportb (port ) ; 
dth=inportb (port+1) ; 
adl=dtl/16; 
adt=dth*16+adl ; 
chv=dtl-adl*16 ; 

/* Reading motor currents, reference speed and rotor speed */ 

switch (chv) 

{case 0: 

ain0=(adt-2048) *5.0/2048; 

break; 

case 1: 

ainl= (adt-2048) *5.0/2048; 

break; 

case 2 : 

ain2= (adt-2048) *5.0/2048; 

break; 

case 3 : 

ain3=(adt-2048)*5.0/2048; 

break; 

case 4: 

ain4=(adt-2048) *5.0/2048; 
break; 

} 

++pool; 

} while (pool<5 ) ; 

} 
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